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A next-generation liquid-scintillator detector will be able to perform high-statistics measurements 
of the solar neutrino flux. In LENA, solar ^Be neutrinos are expected to cause 1.7x10* electron 
recoil events per day in a fiducial volume of 35 kilotons. Based on this signal, a search for periodic 
modulations on sub-percent level can be conducted, surpassing the sensitivity of current detectors by 
at least a factor of 20. The range of accessible periods reaches from several minutes, corresponding to 
modulations induced by helioseismic g-modes, to tens of years, allowing to study long-term changes 
in solar fusion rates. 



I. INTRODUCTION 



Since the early days of solar neutrino flux measurements, efforts have been made to discover a variation 
of the neutrino rate over time. The most prominent effect is the experimentally confirmed annual 
variation of the flux by 7% that is induced by the eccentricity of the Earth's orbit (e. g. [U H])- But 
a variety of alternative sources of flux modulation is conceivable: The survival probability of solar 
electron neutrinos might be influenced by fluctuations of the solar matter density |3"5j or by traversing 
terrestrial matter 6;JS] before reaching the detector on Earth. Solar neutrino production rates might 
even change in the course of the solar cycle of about 11 years [9j, or might be subject to short-term 
variations correlated to the oscillation of the solar core temperature induced by helioseismic waves [TU] . 

The next-generation neutrino observatory LENA will be a low-background liquid-scintillator exper- 
iment [TTHT5] . In recent years, the potential of this detector technique has been demonstrated by the 
first real-time measurement of ^Be neutrinos in the solar neutrino experiment Borexino |14] . Opposed 
to the comparatively low rate of ~50 counts per day (cpd) in Borexino, about 10'* cpd of ^Be neutrino 
events are expected in LENA due to the considerably larger target mass (Sect.|ll|. Given this large 
statistics, it is clear that LENA will be much more sensitive to temporal flux fluctuations as any of the 
preceding solar neutrino experiments, including the Super-Kamiokande detector [15 . In the following 
analysis based on the Lomb-Scargle method [1^1 [T7] (Sect. Ill), we demonstrate that LENA will be sen- 
sitive to fluctuations on a sub-percent level over a wide range of frequencies: Studies of the sensitivity 
cover modulation periods from several minutes up to tens of years (Sect. IV). The resulting discovery 
potential for a number of plausible modulation sources will be presented in Sect.|V| 



II. '^BE NEUTRINO RATE AND BACKGROUNDS 



The low background environment necessary for a low-energy rare-event search is reflected in the 
design of the future LENA detector: The target volume will encompass about 5xl0^m'^ of organic 
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liquid scintillator contained in an inner cylindrical nylon vessel of 100 m height and 26 m diameter. 
Scintillation light is observed by photomultiplier tubes (PMTs) that are mounted to a scaffolding in 
2 m distance from the nylon cylinder: About 45,000 8"-PMTs equipped with light concentrators would 
be needed to reach the intended photocoverage of 30 %. Additional 2x IC* m'^ of non-scintillating buffer 
liquid will surround the scintillator, providing shielding against the 7 rays emitted by the PMTs and 
the surrounding steel tank of 30 m diameter. Outside the tank, an external water volume will serve 
as an active muon veto at the detector perimeter, while the lid of the tank will be covered either by 
plastic scintillator panels or resistive plate chambers (RPCs) to enhance the spatial resolution of muon 
tracking. Further details are given in I13j. 

Depending on the natural radioactivity content of the detector materials used, the fiducial volume 
for solar neutrino detection will have to be chosen considerably smaller than the full scintillator volume. 
Monte Carlo simulations show that the 7 ray background expected in LENA will be dominated by the 
natural radioactivity of the PMT glass in case of a steel tank. The total rate of 7 rays emitted by ''"K 
and elements of the ^ssy g^j^^ 232rp-|^ decay chains is ~170kBq. However, the self-shielding effect of the 
organic liquid reduces the rate of 7 rays reaching the bulk of the detector considerably. 

Tab.|T] shows the fiducial mass Mfid depending on the radius rfid of the chosen cylindric volume. By 
scaling the ^Be-z/ rate of 49 counts per day and 100 1 observed in Borexino, one can derive the total 
expected rate (3'''^ column in Tab.|l]) as function of Mfid that linearly scales with the available volume. 
However, ^Be-v detection in liquid scintillator relies on the detection of electron recoils [T3]. Part 
of these will lie below the detection threshold of 250 keV, corresponding to ~43 % of the total recoil 
spectrum"'^. Therefore, the signal rate remaining inside the observation window from 250-800 keV must 
be compared to the rate of 7 rays reaching the fiducial volume and creating a signal of the same visible 
energy. In the following, a fiducial mass of ~35kt is chosen, corresponding to ^lO'* ^Be-v recoil events 
and ~300 7 events per day inside the detection window (4*'' and 5"^ column). 



Fiducial Volume 
r-fid [m] A/fld [kt] 


Rate [cpd] 


'Be window [cpd] 
''Be-v BG(7) 


13 43.8 

12 36.6 
11.5 33.2 

11 30.1 
10.5 27.1 

10 24.3 


21.5x10^ 
17.9 xlO** 
16.3x10" 
14.7x10" 
13.2x10" 
11.9x10" 


12.3x10" 1.2x10" 
10.2x10" 6.8x10^ 
9.34x10^ 2.5x10^ 
8.42x10^ 7.8x10° 
7.56x10^ 5x10"^ 
6.82x10^ 3x10"^ 



TABLE I: The radius rad of the fiducial volume in LENA as well as the mass contained (Mad) depends on the 
required shielding against the 7 background produced by the PMTs. Signal and background rates are computed 
for LAB (CisHao, p — 0.86 g/^) [18]. The middle column Usts the total rates of ^Be-u events as function of 
Mfld, while the last columns present the rates of ^Be-v events and external 7-rays expected in the ^Be detection 
window from 250 to 800 keV. 

The experience with Borexino shows that the radiopurity requirements of 10^^^ g of ^38^ and ^'^^Th 
per gram of liquid scintillator (or ^^0.5 decays per day and ton) are achievable, which is a prerequisite 
for the detection of the ^Be signal [19]. As the necessary precautions in the production, transport and 
handling of the organic solvents as well as potent purification procedures have been put to test in both 
the Counting Test Facility (CTF) and Borexino itself, similar or even lower levels of radioactivity can 
be expected in LENA. In all of the following analysis, it is assumed that the radiopurity conditions 
achieved in Borexino will be reproduced for LENA. 

Recent developments in the R&D of organic solvents indicate that the optical transparency of the 
liquid to the scintillation light will surpass former expectations [20] . Some brands of linear alkylbenzene 



^ The exact values depend on the energy resolution of the detector, that defines the smearing of the endpoint of the ^"C 
spectrum. 
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— Fit! x'/NDF = 185/174 
— 'Be: 49±3 cpd/100 tons 
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— "C: 25±1 cpd/100 tons 
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FIG. 1: Solar neutrino recoil spectra from elastic electron scattering as measured by Borexino [19]. The 
data points correspond to the measured spectrum below 2 MeV for 15 kton-days of exposure. Both signal and 
background spectra are included in the fit. ^^"Po has been subtracted using a//3 discrimination. 



(LAB) feature attenuation lengths on the scale of 20 m JJJ, corresponding to a light collection efficiency 
of about 450pe/MeV (photoelectrons per MeV) in LENA. This value is only slightly worse than the 
value achieved in Borexino |19j . This is a significant improvement, as self-absorption of the scintillation 
light in the hquid is a major issue and yields of ^lOOpe/MeV would be expected if the technology of 
Borexino was simply duplicated. Due to the enhanced light collection, the low-energy electron-recoil 
spectrum of LENA might look much the same as that observed in Borexino (Fig.[T]) . Therefore, analysis 
techniques based on a fit of pre-calculated signal and background spectra to the overall sum spectrum 
as they are currently used in Borexino are expected to perform equally well in LENA [Mj [19] . 



III. LOMB-SCARGLE PERIODIGRAMS 



A variety of different techniques has been applied for modulation search in solar neutrino data, reach- 
ing from unbinned maximum likelihood methods fl^ [531 to wavelet approaches sensitive to transient 
modulations [1 . In the following Lomb-Scargle (LS) periodigrams are used pQ [TBI [TT] 112 HI] ) a common 
method to analyze a binned data set for periodical modulations of the type 

N{t) Nq- (1 + A-sm{2T:t/T + ip)) . (1) 

Here, N{t) is the expected event rate at time t, while A^'o represents the mean rate, A indicates the 
relative amplitude of the modulation, T describes its period and (p the phase relative to the start of the 
measurement. The power P at which a modulation can be found in the data is obtained by weighting 
the difference between the number of events N{ti) in every data bin i and the expected mean value Nq 
with cosine and sine functions that oscillate with the investigated period T: 




where n is the number of bins, a the standard deviation of N, and U the time at which the data 
corresponding to bin i was acquired [22] . As the quadratic sums of both cosine and sine are used, the 
result is independent of the modulation phase as long as the modulation period is short in comparison 
to the overall measurement time. The weights Wi allow to compensate differences in the bin width w 
of individual data samples. As the analysis presented in Sect. |IV| will be based on equally wide bins 
(neglecting dead time), Wi — w has been used. Therefore, all weight factors cancel each other. The 
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FIG. 2: Lomb-Scargle periodigrams for a MC data set of 2 years measurement time. Left: White noise spectrum. 
Right: A modulation of 2 % relative amplitude and a period of 0.1 years was included. A corresponding peak 
is visible at the expected position, that is also clearly exceeding the white noise level. 



phase factor tp can be derived from the equation 



tan(44)- £r"'^^"l!i! - (3) 



^/ X]iIi'f"»cos(47r|l^ 

Fig. |2] shows on the left a LS periodigram of a white noise Monte Carlo (MC) data set, while on the 
right a modulation oi A = 2% and T ~ 0.1 yrs is included. The LS power P of a given modulation 
is primarily a function of its amplitude A. Statistical fluctuations of the bin content, A — AN/N, 
will alter both the maximum P generated by white noise and the exact value observed for the actual 
modulation. To assess the significance of a modulation discovery, it is therefore necessary to know the 
statistical fluctuations of both the white noise level and the signal height. 

The analyses conducted in the next chapters are based on a common basic scheme that was adopted 
from |22j : To determine whether a modulation of relative amplitude A and period T is pronounced 
enough to be detected, lO"* MC data sets containing a modulation parameterized according to Eq. ^ 
are simulated, and the LS power P for the signal at period T is computed. Then, a white-noise sample 
of 10*^ Monte Carlo (MC) data sets is generated for comparison. The detection threshold is set above 
99.7% of the white-noise LS power distribution, which corresponds to a 3a discovery potential. If 
all of the signal spectrum is above this detection threshold, the probability for a 3a discovery of the 
modulation corresponds to 100 % (compare Fig.js]). In case there is an overlap and the signal distribution 
extends below the detection threshold, the discovery probability is reduced to the proportion of the 
distribution remaining above threshold. In the following analyses, the sensitivity threshold s in LENA 
is defined as the amplitude for which the probability for a 3a discovery is 90 %. 



IV. RESULTS ON SENSITIVITY 



The results presented here are based on a LS study of the expected signal rate generated by ^Be-iy 
electron recoils in LENA. The investigation will be presented in three regions defined by the modulation 
period T: Starting from values of T of the order of days to years in which LENA features the largest 
sensitivity, very short periods of about 30 min and very slow modulations on the scale of centuries are 
explored. A central issue in all analyses is the presence of radioactive background that superimposes 
the neutrino signal. The strategy to isolate the signal will differ depending on the chosen bin width w 
and will be discussed in the individual sections. 
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FIG. 3: Distribution of Lomb-Scargle powers of white-noise (solid line) and of a modulation of 2 % relative 
amplitude and 0.1 years period (shaded area) (compare Fig.[2|. There is a clear separation between the distri- 
butions. The sensitivity threshold s is set at 99.7% of the white noise distribution (dashed line). 



A. Medium- Term Modulations 



First, the dependence of the detection sensitivity on measurement time and bin-content uncertain- 
ties will be outlined. Then the uncertainty level that can be expected from the ^Be-v analysis will 
be discussed, relying on the experience gathered in Borexino. Finally, the resolution of neighboring 
modulation lines in the LS periodigram will be discussed. 



Primary Parameters Governing Sensitivity 



The Lomb-Scargle method is most sensitive for modulations in an intermediate range of periods: T 
should not significantly exceed the overall measurement time r, but should also not be so short to come 
into conflict with the bin width w. MC tests show that the sensitivity is independent of T and also of 
the modulation phase tp if these conditions are fulfilled (see also Fig.[6|. This aspect will be discussed 
in Sect. UVCl 

One of the main factors determining the minimum detectable amplitude A is the measurement time r. 
Moreover, the relative uncertainty of the event number in each bin, A = AN/N, has a large influence. 
To demonstrate these dependences, Tab.[II| shows LENA's sensitivity as a function of t and A. The 
chosen bin width w — 1 day. As stated in Sec. Ill s corresponds to the minimum amplitude that is 
discovered at 3(t level in 90 % of all MC data sets. Samples of 10 000 MC data sets each were used to 
determine the detection threshold by comparison of white noise and signal spectra. 



Uncertainties from the Fit 



To assess the expected sensitivity s in LENA, it is necessary to determine both statistic and systematic 
uncertainties of the neutrino-rate measurement during a period corresponding to the bin width w. As 
the ^Be-z/ rate surpasses all other neutrino rates by at least one order of magnitude, these neutrinos 
are the most sensitive probes for this analysis ^. Based on a fiducial mass of 35 kt, the expected ''Be 



^ The question whether neutrinos might be another promising candidate for measurements of solar temperature 
fluctuations will be discussed in Sect.fVl 
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[yrs] 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 



Relative uncertainty per bin A = AN/N [%] 
123456789 10 



0.34 
0.26 
0.19 
0.19 
0.18 
0.17 
0.17 
0.15 
0.14 
0.11 



66 0, 
47 0, 
38 0, 
34 0, 
29 0, 
28 0, 
28 0, 

25 0, 

26 0, 
21 0, 



97 1.29 
68 0.91 
57 0.75 
49 0.65 
45 0.58 
40 0.53 
39 0.49 
36 0.49 
36 0.44 
31 0.41 



1.61 1, 
1.14 1, 
0.93 1, 
0.80 0, 
0.72 0, 
0.66 0, 
0.61 0, 
0.57 0, 
0.54 0, 
0.50 



93 
36 
19 
97 



2.25 
1.59 
1.29 
1.12 
86 1.08 
79 0.92 
73 0.85 
69 0.80 
65 0.80 
62 0.71 



2.57 

1.80 

1.48 

1.28 

1.15 

1.05 

0.98 

0.91 

0.92 

0.81 



92 3.21 
03 2.26 
66 1.84 
45 1.60 
29 1.43 
18 1.31 
09 1.20 
1 
1 



02 
97 



.14 
.15 



91 1.02 



TABLE II; Sensitivity s of modulation search as a function of measurement time r and of the relative uncertainty 
of the rate measurement A = AN/N. 10* MC data sets were generated. The sensitivity s corresponds to the 
minimum amplitude (in %) that would be observed at 3a" level for 90 % of the MC data sets. 



neutrino rate will be of the order of 10'* cpd. Consequently, a statistical error of 1 % can be expected 
for a standard bin width w = Id. 

However, the data collected in Borexino suggests that there are several sources of radioactive back- 
ground inside the ''Be detection window. The most prominent ones are the decays of ^^°Po, ^*°Bi, and 
^^Kr dissolved in the scintillator. 

The isotope ^^"Po emits a-particles with an energy of 5.41 MeV. Due to the quenching effects present 
in liquid scintillators, this line is shifted to about 450 keV in the electron recoil-spectrum (depending on 
the specific scintillator composition) [TH |T5] . Assuming for LENA a replication of the initial radioactive 
background levels of Borexino, the contamination with ^^''Po would be considerable, in the beginning 
surpassing the expected ^Be rate by more than two orders of magnitude. Despite this large rate, 
the effect of ^^°Po on the ""Be rate determination can be assumed as negligible: As demonstrated in 
Borexino, the efficiency of liquid-scintillator detectors to identifiy a events by a pulse-shape analysis 
reaches levels of at least 98% Moreover, the ^^"^Po half-life of 138 d is short compared to the 

expected operation time of LENA. The initial value of ^*°Po contamination will be reduced each year 
by more than a factor of 6 as long as it is not fed by the decay of long-lived ^^°Pb. 

The ultimate limit for the ^^°Po reduction by decay is set by the ^^°Pb contamination of the scin- 
tillator, which produces ^^°Po by the beta decay chain ^*°Pb ^^"Bi — ^ ^*°Po. In Borexino, the 
initial amount of ^^"Po by far surpassed the concentrations of ^^'^Bi and ^^"Pb; the ^^^Rn chain is out 
of equilibrium. However, ^^°Pb with a half-life of 22.3 yrs is relatively long-lived. Even assuming an 
operation time of 30 years for LENA, the final ^^'^Pb will still be at about 40 % of the start value. 
While the /3 spectrum of the ^^"^Pb decay itself features a sufficiently low endpoint of 64keV not to 
be in conflict with neutrino detection, its comparatively short-lived decay daughter ^*°Bi features a /3 
endpoint of 1162 keV, severely interfering with the signals of ^Be, pep, and CNO neutrinos. As the 
emitted particle is an electron it cannot be distinguished from a neutrino-induced electron recoil by its 
pulse shape. The only possibility is an identification of this contribution to the overall recoil spectrum 
by its known spectral shape. Similar arguments can be made in the case of ®^Kr, a /3-unstable isotope 
with an endpoint energy of 687 keV and a half-life of 10.8 yrs. 

For an estimate of the impact of these backgrounds and other systematic sources on the overall 
uncertainty in LENA, we refer to the results of Borexino. The ^Be neutrino analysis presented in [19] is 
based on a data set corresponding to 198 live days of data taking or about 7500 neutrino events. This 
corresponds to 1/2 of the statistics collected in LENA in a single day. As described before, the ''Be 
neutrino rate for the overall period was determined by a spectral fit to the data, taking into account 
the recoil spectrum of ^Be and the contribution of ^^Kr and ^^°Bi background events. For Borexino, 
the resulting ''Be event rate is 49 ± 3(stat) i 4(-syst) cpd. The uncertainties are mostly caused by the 
determination of the fiducial mass and by the detector response function. 

Unlike the uncertainty of the absolute ^Be neutrino rate, a relative rate measurement as in the case 
of the modulation search will not be affected by the exact knowledge of the fiducial mass. On the other 
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hand, the uncertainty in the detector response function enters the analysis by ambiguities in the spectral 
fit. It can be significantly reduced by calibrations based on the insertion of radioactive sources in the 
detector. Both Borexino and Kamland have undertaken dedicated campaigns. Based on the result of 
the calibrations, the Borexino collaboration is now aiming at a reduction of the uncertainties on the 
^Be measurement, corresponding to an uncertainty of ^3 % for the ''Be rate extracted from spectral fits 
and an overall uncertainty of ~5 % including other systematics. In principle, the uncertainty obtained 
in LENA for a single-day measurement might be even lower than this value: The average background 
rate per day can be determined using the statistics corresponding to time spans of months or years, and 
the remaining uncertainty on the fiducial mass can be neglected. Excluding a change in background 
rates due to novel contaminations during the measurement time, the level of radioactive contamination 
of the scintillator is constant (ignoring for the moment the decay of radioimpurities) . Using the mean 
background rates as pull parameters for the spectral fit will further reduce the uncertainty of the 
day-by-day determination of the ^Be rate. 

Conservatively, an uncertainty of A = AN/N = 3 % can be expected for the individual bins, corre- 
sponding to 10'* events as the bin width is a day. Assuming a measurement duration and a stability 
of the fluctuation of 1 (lO)yr(s), a sensitivity of about 1%(0.3%) can be reached for medium-term 
fluctuations (compare Tab.[ll|. 



Resolution of Modulation Period 



It is possible that the solar neutrino flux is subject to multiple superimposed modulations, featuring 
either a whole spectrum of frequencies or several discrete values. In the latter case, the ability to distin- 
guish between nearby lines in the LS periodigram and the determination of the statistical uncertainty 
in the peak frequency are of great interest: 

The minimum difference in modulation period that is needed for resolving two close-by superimposed 
modulations is determined by the corresponding line width AT in the LS periodigram. In Fig. [2] the 
relative width AT/T is 2 % (Icr) for a modulation of T = 0.1 yrs and t — 2 yrs. 

MC calculations show that both the measurement uncertainty per bin A and the modulation am- 
plitude A have no influence, at least as long as A > s (Tab. [III). Merely the measurement time r and 
the period T alter the resulting AT. Varying both parameters, AT proves to be inversely proportional 
to T. At the same time, it scales as the square of the period T^. This behavior is equivalent to the 
assumption that only the number of available measurement bins n = T/w — T /[t /T) = T"^ /t has an 
influence on the accuracy. Assuming a bin width of 1 day, AT/T can therefore be written as: 



T yO.lyr/ \ t J Vly/ V''" 

In case of a single modulation, the mean value of the line, i. e. the modulation period T, can be de- 
termined at high accuracy from the LS periodigram (Fig.|2]). To determine the jitter in the maximum 
value of the line, periodigrams of 100 MC data sets were created according to the parameters mentioned 
above. The peak value featured a root mean square of 6xl0~^ relative to T. 



B. Long-term Changes in Flux 

The LS method is most sensitive to the discovery of flux oscillations for periods shorter than the 
measurement time. The currently planned operation time of the LENA detector is 10 to 30 yrs, which 
is of the order of magnitude of a full solar cycle. However, it is not excluded that there might be 
modulations of the total nuclear energy production of the Sun on longer time spans which would be 
reflected in the emitted neutrino flux. Therefore, a series of MC calculations was performed to determine 
the sensitivity to very long-term modulations dependent on the period T. 

Different to the situation for medium-term modulations, the sensitivity to slow oscillations also is 
a function of the modulation phase ip relative to the start of the measurement period. For t <^ T, 
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FIG. 4: Influence of the modulation phase ip on the sensitivity s to a modulation of T = 200 yrs, that varies by 
an order of magnitude (t = 10 yrs, A = 3%). 




FIG. 5: Influence of the modulation phase i/p on the sensitivity s for modulations with periods longer than the 
measurement time, T > t {t — lOyrs, A = 3%). 



the visible change in detected neutrino rate will be smallest near the flux maxima and minima, N = 
Nq{1 ± A), while it will be largest near the zero points of the oscillation function, N ^ Nq. This 
expectation is confirmed by MC calculations displayed in Fig.|4] For a measurement time t — 10 yrs 
and A — AN/N = 3%, the minimum detectable amplitude A oi a modulation with a period of 200 yrs 
varies between 1.3% and more than 17%, depending on ip. 

The predicted sensitivity depends on the actual combination of r, T, and ip. This is illustrated in 
Fig. [5] that for r — 10 yrs scans the range of modulation periods from 10 to 500 yrs for various values 
of if. For T > 2r, the phase (p becomes the dominant factor in determining the achievable sensitivity, 
while for T > 5r a disadvantageous value of (p is able to obscure a modulation that would otherwise be 
well visible. The longer the operation time of the detector, the more these effects will be reduced. 



C. High- Frequency Modulations 



While the main concern for LS modulation search at low frequencies is the total measurement time r, 
the central issue for very high frequencies is the bin width w and the corresponding uncertainty on the 
bin contents A{w). As discussed in Sect. Ill there are certain techniques as e. g. unbinned maximum- 
likelihood methods that do not suffer from the the finite bin width [22]. But also in this case, statistical 
subtraction of the background will limit the achievable sensitivity. 
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FIG. 6: Sensitivity dependence on modulation period T versus number of bins n. The plot is generated assuming 
n — 100 and A = 1 %. The minimum period Tmin that still provides an optimum sensitivity corresponds to the 
width of five bins. 



Minimum Length of Modulation Period 



The LS method loses in sensitivity s in the case that the modulation period T gets short compared 
to the bin width w. Fig. [6] shows s for a 100 bin MC data set with a bin content uncertainty A = 1 %. 
As mentioned already in Sect. IV A the sensitivity is almost constant over a wide range of periods, 
5w < T < 150it;. The minimum period detectable at optimum s is therefore 



5 w. 



(5) 



Statistical Uncertainties for Short Time Binnings 



As shown in Sect. IV A the ^Be-iy rate in LENA is best determined by a spectral fit for bin widths 
w > 1 day. However, this method will fail for w I day as at some point the fit will return less precise 
results than a simple counting of the events inside the bin. If this counting approach is chosen, the 
background contributions inside the bin must be statistically subtracted from the total bin content A'tot 
in order to obtain the signal rate. As discussed in Sect. |IV A} the primary background contributions 
are ^^Kr and ^^"Bi; ^^°Po can be removed by pulse-shape discrimination (or decay). The number of 
^Be-j/ events in a single bin, iV^(u)), is therefore given by: 

N,iw) ^ Ntotiw) - Nkg{w), (6) 

where A'^^g — A'xr + A^Bi- The number of background events scales with the background rate rbg and the 
bin width w, N^g = r\^gW. The total number of events is A'tot = (j'i^+r]^g)w. In the following assessment, 
background rates are taken corresponding to the situation in Borexino [14]: The signal-to-background 
ratio inside the detection window is roughly 1:1, so the rates in LENA are 

« rbg « lOVd. (7) 

The contribution of background events to A^tot varies according to the statistical uncertainty y^A^bg- 
Moreover, the statistical subtraction of A'bg from TVtot is based on rbg. The mean rate is determined 
by a spectral fit applied to a larger data sample, collected over a time span in which the contamination 
remains unchanged. Very conservatively assuming that the value of rbg is determined on the basis 
of 1-day data sets, the corresponding uncertainty is of the order of 3 % (as it is the case for ^Be-z^'s 



of Sect. IV A I. Therefore, the total uncertainty on the neutrino rate per bin, AN^^w), is determined 



10 




bin length (days) 



FIG. 7: The relative uncertainty of the bin content of ^Be events depending on the bin width w. The 
numbers are based on the statistical subtraction of the background rates caused by *^Kr and ^^"Bi decays. 



by the quadratic sum of the uncertainty for the total number of events, AiVtot — \J ^tot , and for the 
background contribution AA'bg — \/ N^g + (O.OSA'bg)^: 

AN, = ■^(AArtot)2 + (AiVbg)2 = ^Ntot + A^bg + (O.OSTVbg)^ ^^{r, + 2r^^)w + (O.Oarbg)^^^. (g) 

By entering the rates quoted in ([7|, the relative uncertainty of the ^Be-i/ rate per bin can be written 
as: 



AN, AN, , /2.9d 

N, r,w V w 

For w = Id, A, = AN,/N, « 3% is governed by the uncertainty in the determination of the 
background rate. But for bin width w < 0.3 d, the statistical uncertainties of the individual bins 
will dominate the uncertainties of the long-term background rate determinations. Fig. [7] shows the 
corresponding function of the uncertainty, A,{w). 



Single-Bin Uncertainty vs Number of Bins 

A fixed relation between the sensitivity s and the minimum measurement time r can be obtained by 
combining the bin content uncertainties A and the number of bins n for the LS analysis. In accordance 
to Eq.[5j only sufficiently long periods Tmin > 5 if are considered. Fig. [8] shows the results of MC 
simulations describing the dependence of s on n. As one would expect, the sensitivity of the analysis 
improves with n. In addition, the uncertainty of the bin content A was varied. From the MC result, it 
is possible to derive a function describing the dependence of s on n and A: 

s{n,A) ^m - A-n"^. (10) 

The actual fit parameters vary for the individual values of A and are listed in Tab. |III[ Approximately, 
the sensitivity can be described by 



s(n,A)sa7.5-A-n~'' 



(11) 
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FIG. 8: Dependence of the sensitivity s (Tmin > 5m) on the number of bins n available for LS analysis. The 
markers represent the results obtained by simulation, varying both the bin-content uncertainty A and n in the 
MC data sets. The lines correspond to fits to the MC results at distinctive values of A, following Eq. ( |10[ ). 
Parameters are listed in Tab.lIIII 



A 


m 


7 


0.03 


7.4847(13) 


0.51880(3) 


0.05 


7.3248(4) 


0.516292(13) 


0.1 


7.5843(8) 


0.520261(19) 


0.15 


7.550(2) 


0.52059(5) 


0.2 


7.4675(15) 


0.51792(4) 


0.3 


8.605(2) 


0.5403(5) 


0.4 


7.875(7) 


0.52803(18) 


0.5 


7.938(6) 


0.53077(14) 



TABLE III: Fit parameters for the sensitivity s as function of the relative bin content uncertainty A and the 



number of measurement bins n, following Eq. ( 10 1 



Sensitivity Estimate 



The relation shown in Eq. ( 10 ) describes the dependence of the sensitivity s of a modulation search 
on the bin content uncertainty A and the number of bins n. For estimating the final sensitivity of 
LENA, the interesting parameters are the minimum modulation period T^^n detectable as a function 
of the overall measurement time r. 

As stated in Eq. ([s]), Tmin has to include at least 5 bins. The relation between A and the chosen bin 
width w is described by Eq. ([9| . The measurement time r is the product of bin width and number of 
bins, T = w ■ n. The result of combining this information with Eq. (10 1 is depicted in Fig.[9j The best 



sensitivity s achievable is s « 2 % for a measurement time of 1 month, and s « 0.5 % for r = 1 yr. The 
worsening of s with increasing T visible in Fig. [o] is based on the shape of Eq. ([9|: For large T, AA^ 
is dominated by the constant term originating from the long-term background fit. However, it can be 
avoided by choosing smaller bin widths w <^ T even in the case of larger modulation periods, as Eq. ([s]) 
only provides a lower limit. Therefore, the minimum sensitivities quoted above are constant as long as 
T < T (compare Fig. [6]). 

The sensitivity levels obtained by this method would surpass the values presented in Sect. |IV"A| 
provided sufficient measurement time can be reached. If one considers a data set of n = 3.65 • 10^ bins 
of ui = 10^^ d in width (corresponding to 10 years of measurement), the corresponding bin uncertainty 
is A^ « 17.3% from Eq. ([9]). In this case, Eq. (11) returns a sensitivity s « 0.2%, close to the 



maximum sensitivity s — 0.1% obtained if systematical errors are not considered (corresponding to 
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FIG. 9: Sensitivity s in LENA for shor t m odulation periods. The plot is generated using the dependence of 
sensitivity on the number of bins of Eq. ( 10 1. Each data point corresponds to an individual bin uncertainty and 
the correlated minimum modulation period that can be detected. For small periods, the greater number of bins 
is leveled by the increase in bin uncertainty and the correlated loss in sensitivity. 



A = 1 % in Tab.|TT]). However, the presented analysis relies on the uniformity of both detector response 
and background levels over the whole measurement period. The experience of Borexino shows that 
this will be hard to achieve on the scale of several years. Shorter stability periods can be much easier 
accommodated in an analysis following Sect. |IV A} 



V. APPLICATIONS AND DISCUSSION 



Based on the sensitivity estimates presented in Sect. |IV[ one can define three modulation period 
ranges that will be accessible in LENA: 

• 10 min < T < Id: For very short periods, the LS analysis relies on simple counting of events 
inside the neutrino detection window. The sensitivity mainly depends on the measurement time 
r: For r = 1 yr, s sa 0.5 % in amplitude can be reached. 

• Id < T < lOyrs: For medium-range periods, s « 0.3% can be achieved for a measurement time 
of lOyrs based on fits to the electron recoil spectrum in day-long bins. 

• T > lOyrs: For long periods, the sensitivity decreases with T, and is also subject to the strong 
influence of the phase shift ip relative to the acquisition time window. Periods beyond a few 
hundred years are generally not accessible. 

In the following, several possible sources of modulations in the solar neutrino flux and the sensitivity 
of LENA for their detection will be discussed. 



Keplerian Motion 



Although already discovered in preceding solar neutrino experiments, an analysis would return a 
confirmation of the annual flux modulation caused by the ellipticity of the terrestrial orbit at high 
significance (compare Tab.[Tl|). For short measurement times, t « lyr, it could be argued that the 
presence of the large signal might cause a decrease in sensitivity for neighboring frequencies in the LS 
periodigram. However, MC simulations investigating this effect show that an artificial suppression of 
this signal in a data set by a bin-wise division with A{t) 0.04sin(27r(t/l yr)) (see Eq. ([T])) compensates 
the distortions in the baseline of the periodigram. 



13 



Day-Night Effect 

The observation of a day-night asymmetry in the solar neutrino fiux, induced by the transition of 
solar neutrinos through the Earth's matter, has been considered as a "smoking gun" evidence for the 
MSW effect [6 . Formerly, the LOW solution for neutrino oscillation parameters predicted rather large 
amplitudes of up to 40% for the observable ''Be flux [71[S]. Due to the consolidation of the MSW-LMA 
model of neutrino oscillations, such a large effect seems now very unlikely. Alternatively, a similar 
day/night asymmetry of the survival probability is predicted in the case of mass varying neutrinos: the 
expected amplitude is ~20%, while the phase is opposed to regular matter effects [55]. Up to now, 
no indication of modulations at this scale has been observed Borexino data. Nevertheless, it might be 
worthwhile to search for such a modulation in LENA "^Be data. Due to the long-term stability of the 
effect, sensitivities of 10"'^ could be reached for a diurnal modulation. 

Correlations to the Solar Cycle 

A variety of authors has analyzed the data collected over long measurement periods in the radiochem- 
ical experiments, searching for low-frequency modulations with periods of the order of tens of years 
(e.g. [27J). In the data reported by all gallium experiments (Gallex/GNO and SAGE), a decrease in the 
neutrino flux on the order of 20 % over about 10 years of combined measurement time is visible [24ll28j. 
However, the statistical uncertainties on the individual measurement bins are too low to provide a 
significant result; data is as well consistent with a constant neutrino rate. 

In terms of the analysis presented above, the inconclusive change in flux could be approximated by 
a modulation of T > 20 yrs and A= Confirmation or rejection by LENA is therefore well within 

reach, even if the actual measurement period were only several years. However, an evaluation of the 
periodicity of the signal would probably require a measurement time of the order of 10 years. 

Helioseismic Waves in the Neutrino sphere 

It has been argued that helioseismic waves might give rise to modulations in the solar neutrino 
production rate jlOj . These waves are present throughout the Sun and stretch over a large range of 
amplitudes and frequencies |29| . Pressure-driven waves (p-modes) propagate close to the solar surface 
and have been extensively studied by optical helioseismology. However, buoyancy or gravity-driven 
waves (g-modes) are confined to the inner regions of the Sun due to their dampening in the convective 
zone. Up to now, there is no unambiguous evidence for g-modes by standard helioseismology, although 
hints of an oscillation at a frequency of 220.7 jrHz have been found in SOHO data [3D]. 

The expected frequency range of solar g-modes can be derived from solar model calculations: Typical 
values are 200]jHz or lower, corresponding to periods longer than 1.4 h. The current best limit for a 
corresponding modulation in the observed neutrino flux has been derived from SNO data and excludes 
amplitudes of more than 10% [lOj. Based on the present study, LENA will be able to surpass the 
accuracy of this measurement by at least a factor of 20 in the relevant frequency region. 

There is no specific mechanism that relates g-mode amplitudes to a modulation of the neutrino fiux. 
However, the primary possibilities are a change of neutrino production rate induced by a localized 
change in temperature, or a variation of the matter potential infiuencing the MSW survival probability 
(see below). 

In the first case, the sensitivity of the search would be enhanced by the steep dependence of the 
fusion rate on the temperature of the medium |31) . Standard solar model calculations indicate that the 
^Be rate scales with the 11th power of solar core temperature. This implies that relative temperature 
changes on the level of s/11 w 5 x 10""* might be detectable in LENA. The variation will be even 
more pronounced in the *B fusion rate, as the temperature dependence <j)»^ oc T^^^ is even steeper. 
The considerably lower sensitivity (compared to that of ''Be-i') to ^B-j/ fiux modulations of s ~ 1 % in 
LENA, corresponding to a detection rate of ^100 cpd, would be roughly compensated by the enhanced 
temperature sensitivity. 
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Helioseismic waves are usually quantified by means of spherical harmonics. Solar model calculations 
indicate the largest oscillation amplitude for g-modes of low quantum numbers [29j . This is fortunate, 
as the full effect of temperature variation will only be observed if all of the ^Be (or ^B) fusion region 
undergoes the oscillation in unison. 

Solar Matter Effects and Spin Conversion 

The influence of the solar matter on the survival probability of *B neutrinos is a basic ingredient 
of the MSW-LMA oscillation model. However, the predictions differ concerning the question whether 
changes in solar density, as induced by helioseismic activity, are sufficiently pronounced to influence 
the survival probabilities [32*. 

While LENA'S primary sensitivity for modulation search is in the ''Be signal, it seems very unlikely 
that matter effects will have a significant influence on the survival probability at sub-MeV energies. 
However, the ®B neutrino rate will be sensitive to such effects. For a rate of 100 cpd, a minimum 
amplitude of 1 % could be identified. This surpasses the current best limit held by SNO by an order of 
magnitude 10]. 

Another possible cause for modulations in the solar i/g flux is the conversion of neutrino flavors, 
or even neutrino-antineutrino conversion, based on the interaction of a non-zero magnetic moment of 
neutrinos with the solar magnetic fleld |33| . The flavor conversion was originally an alternative oscilla- 
tion mechanism to mass-induced oscillations and is disfavored as the dominant effect by experimental 
results. However, an analysis monitoring the survival probability as a function of the solar magnetic 
field strength is likely to return similar sensitivities as the periodic modulation searches. It will also be 
sensitive to the v ^ v conversion. 

For the latter case, both KamLAND and Borexino have performed direct searches for the appearance 
of Df, events in the energy region above 8 MeV, in which the overwhelming background of reactor Vf. 
declines [34, 35 . It is obvious that LENA could put much more stringent limits due to its large target 
mass. However, also in this energy region there is an irreducible background due to the component 
of the diffuse Supernova neutrino background, which poses a natural limit |36j . 

VI. CONCLUSIONS 

The next generation of large-volume, low-energy neutrino detectors will enhance the experimental 
sensitivity towards modulations in the solar neutrino flux [37j . A liquid-scintillator detector seems 
especially attractive due to the possibility to observe the ''Be neutrino line at 866 keV: As the expected 
event rate is of the order of lO'' counts per day, the high statistics will allow to search for modulations on 
a sub-percent level, by far surpassing the sensitivity of currently running experiments. The sensitivity 
of the search stretches from very short time scales of the order of tens of minutes to tens or even 
hundreds of years. This will allow to probe the frequency regions of the helioseismic observation for 
g-mode oscillations in the solar center [53] , the day- night effect [Tj H] , and to search for variations of 
the fusion rate with the solar cycle. 
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